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GUIDANCE AND CONTROL STRATEGIES FOR AEROSPACE VEHICLES 

(NASA Langley Grant NAG1-736) 


PART I 


This final report under NASA Langley grant NAG1 736 consists of two 
parts. Part I consists of the summary of the earlier work whose reports have 
been submitted at various times during the period (1/1/87 to 12/31/1989). Part 
II describes in detail the research work done during the period 1/1/1990 to 
7/31/1990. 
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I. INTRODUCTION 


The first part of the report concerns broadly the summary of the work 
done in the areas of singular perturbations and time scales (SPaTS), 
aerobraking technology, guidance and aerocruise. 

(i) SPaTS 

The dynamics of many control systems is described by high-order 
differential equations containing parameters such as small time constants, 
masses, moments of inertia, inductances, and capacitances. The presence of 
these "parasitic" parameters is often the source for the increased order and 
the "stiffness" of the system. The "curse" of the dimensionality coupled with 
the stiffness poses formidable computational complexities for the analysis and 
control of such large systems. Singularly perturbed systems are those whose 
order is reduced when the parasitic parameter is neglected. The methodology of 
singular perturbat ions and time scales (SPaTS) is a "boon 1 to control 
engineers in tackling these large scale systems. As such it is very desirable 
to formulate many control problems to fit into the framework of the 
mathematical theory of SPaTS. 

The methodology of SPaTS has an impressive record of applications in a 
wide spectrum of fields including flight mechanics and trajectory 
optimization. The aerospace problems involve, in general, the solution of 
nonlinear differential equations by resorting to numerical integration. 
Analytical solutions are important in providing a general understanding of the 
structure of solutions and a better foundation for the solution of guidance 
problems. With this in view, attempts have been made to obtain approximate 
analytical solutions for the atmospheric entry problem using asymptotic 
methods such as the method of matched asymptotic expansions, singular 
perturbation method, and multiple scale method. 

In this report, using the theory of SPaTS, the various types of 
aerospace-related problems investigated were digital flight control systems 
[1-3, 6-8, 10, 14, 16] and atmospheric entry problems [9-11, 18, 21, 22]. 


* The numbers in brackets indicate the item under List of Publications. 
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(ii) Aerobraking Technology 


The specification spectrum for the proposed Space Transportat ion System 
(STS) places heavy emphasis on the development of reusable avionics subsystems 
having special features such as vehicle evaluation and reduction of ground 
support for mission planning, contingency response and verification and 
validation. According to the report of the National Commission on Space, 
PIONEERING THE SPACE FRONTIER, the concept of aerobraking for orbit transfer 
has been recognized as one of the critical technologies and recommended for 
demonstration projects in building the necessary technology base for 
pioneering the space frontier. In space transportation systems, the 
aerobraking, defined as the deceleration resulting from the effects of 
atmospheric drag upon a vehicle during orbital operations, opens new mission 
opportunities, especially with regard to the initiation of a permanent space 
station. 

The main function of space transportation system is to deliver payloads 
from Earth to various locations in space. Until now, this function has been 
performed by various rockets, the space shuttle, and expendable upper stages 
using solid or liquid propellants. In particular, considering the economic 
benefits and reusability, an orbital transfer vehicle (OTV) is proposed for 
transporting payloads between low Earth orbit (LEO) and high Earth orbit 
(HEO). The two basic operating modes contemplated for OTV are a ground-based 
OTV which returns to Earth after each mission and a space-based OTV which 
operates out of an orbiting hanger located at the proposed Space Station 
Freedom. 

In a typical mission, a space-based OTV, which is initially at the space 
station orbit (SSO), is required to transfer a payload to geosynchronous Earth 
orbit (GEO), pick up another payload, say a faulty satellite, and return to 
mate with the orbiting hanger at SSO for refurbishment and redeployment of the 
payload. The OTV on its return journey from GEO to SSO needs to dissipate some 
of its orbital energy. This can be accomplished by using an entirely 
propulsive (Hohmann) transfer in space only or a combination of propulsive 
transfer in space and aerobraking maneuver in the atmosphere. It has been 
established that a significant fuel savings and hence increased payload 
capabilities can be achieved with propulsive and aerobraking (or aeroassisted) 
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maneuvers instead of all-propulsive maneuvers. The word "aeroassisted" is 
often used to convey that the atmosphere is used to achieve the desired 
deceleration. This leads to an aeroassisted orbital transfer vehicle (AOTV), 
which on its return leg of the mission, dips into the Earth s atmosphere, 
utilizes atmospheric drag to reduce the orbital velocity and employs lift and 
bank angle modulations to achieve a desired orbital inclination. Basically, 
the AOTV performs a synergetic maneuver, employing a hybrid combination of 
propulsive maneuver in space and aerodynamic maneuver in the atmosphere. 
Broadly speaking, the two kinds of orbital transfer are coplanar orbital 
transfer and noncoplanar orbital transfer (or orbital transfer with plane 
change) . 

In this report, using algorithms based on industry standard program to 
optimize simulated trajectories (POST), and multiple shooting method, 
investigations have been carried out to generate fuel-optimal trajectories for 
coplanar orbital transfer [4, 12. 17, 25], and noncoplanar orbital transfer 

[13, 15, 19, 20, 29] arising in aerobraking technology. 

(iii) Guidance 

An optimal trajectory is computed for a given nonlinear dynamical system 
with a fixed set of conditions. However, variation of the initial and final 
conditions, plant parameters would alter the optimal trajectory. It is 
computationally tedious and expensive to repeat the whole optimization 
procedure for every changed condition and obtain a new optimal trajectory. In 
such a situation, an alternative is to linearize the original system and 
generate an optimal trajectory in the neighborhood of the original optimal 
trajectory, involving considerably less computational effort. 

Guidance is the determination of a strategy for following a nominal 
flight in the presence of off-nominal conditions, wind disturbances, and 
navigation uncertainties. In a typical guidance scheme, the final steering 
command is generated as the sum of two components, an open-loop actuating 
(control) signal yielding the desired vehicle trajectory in the absence of 
external disturbances, and a linear feedback regulating signal which reduces 
the system sensitivity to unwanted influences on the vehicle. 
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In this report, guidance schemes for atmospheric maneuver for both 
deterministic and stochastic cases have been investigated [23, 24, 27]. 

(iv) Aerocruise 

There are basically three methods of plane change, (i) impulsive method, 
(ii) aeroglide method, and (iii) aerocruise method. In impulsive method, the 
plane change is achieved entirely outside the atmosphere, and fuel consumption 
is prohibitively large for sizable changes of orbital plane. In both aeroglide 
and aerocruise methods, rockets are used to deflect the vehicle into the 
atmosphere, and the plane change is accomplished by heading change of the 
vehicle. With aeroglide there is no thrusting during the atmosphere, and with 
aerocruise, atmospheric drag is balanced by a continuous thrust to keep the 
spacecraft at a constant altitude and velocity. Propellent expenditure 
comparisons among the three methods of plane change show that the aerocruise 
method is superior to other competing methods for plane changes greater than 
about 20 degrees, and with heating restraints. The basic effect of propulsion 
during aerocruise is to (i) balance drag in order to maintain constant 
velocity, (ii) augment lift with a component of thrust, thus increasing 
cruising altitude over what it would be during aeroglide turn, and finally 
(iii) provide a lateral component of thrust giving the required turn necessary 

for plane change. 

In this report, research has been conducted into cruise maneuver being 
performed using either bank control with constant thrust, or thrust control 
with constant bank control [26, 28]. 
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II. SUMMARY/ABSTRACT OF RESEARCH WORK 

(NASA Langley Grant NAG1-736) 


7 


SINGULAR PERTURBATION AND TIME SCALE 

* 

APPROACHES IN DISCRETE CONTROL SYSTEMS 


Dr. D. S. Naldu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 

and 

Dr. D. B. Price 
Spacecraft Control Branch 
NASA Langley Research Center 
Hampton, VA, 23665 


Abstract: The theory of singular perturbations and time scales (SPaTS) has 
been a powerful analytical tool in the analysis and synthesis of continuous 
and discrete control systems. In this paper, we first consider a singularly 
perturbed discrete control system. Using singular perturbation approach, outer 
and correction subsystems are obtained. Next, by the application of time scale 
approach via block diagonal lzation transformations, the original system is 
decoupled Into slow and fast subsystems. To a zeroth order approximation, the 
singular perturbation and time scale approaches yield equivalent results. 
Roughly speaking, the zeroth order approximation Is somet imes called the 
first approximation. This result is similar to a corresponding result in 
continuous control systems. 


* See items 1 and 24 under List of Publications 
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SINGULAR PERTURBATIONS AND TIME SCALES 

* 

IN DISCRETE CONTROL SYSTEMS- AN OVERVIEW 


Dr. D. S. Naldu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 


Dr. D. B. Price 
Spacecraft Control Branch 
NASA Langley Research Center 
Hampton, VA, 23665 

and 

Dr. J. L. HI bey 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 


Abstract: This paper presents an overview of recent developments in the theory 
of singular perturbations and time scales (SPaTS) in discrete control systems. 
The focus is in three directions: modeling, analysis, and control. First, 
sources of discrete models and the effect of the discretizing interval on the 
modeling are reviewed. Then the analysis of two-time scale systems is 
presented to bring out typical character 1st lc features of SPaTS. Finally, in 
the control of the two-time scale systems, we the Important issue of 
multirate sampling is addressed. The bibliography containing over 100 titles 
is included. 


* See items 2 and 8 under List of Publications 
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# 

ON THE METHOD OF MATCHED ASYMPTOTIC EXPANSIONS 


Dr. D. S. Naldu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 


and 


Dr. D. B. Price 
Spacecraft Control Branch 
NASA Langley Research Center 
Hampton, VA, 23665-5225 


Abstract: A critical examination of the method of matched asymptotic 
expansions (MAE) reveals that the various terms of the common solution ol MAI. 
can be generated as polynomials in stretched variable without actually solving 
for them from the outer solution as is done presently. This also shows that 
the common solution of the method of MAE and the intermediate solution of 
singular perturbation method are the same and hence that these methods give 
identical results for a certain class of problems. An illustrative example is 

given. 


* See items 3, 7 and 18 under List of Publications 
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IMPACT OF ATMOSPHERIC DENSITY SCALE HEIGHT ON THE PERFORMANCE 
OF AEROASSISTED COPLANAR ORBITAL TRANSFER VEHICLES 


Dr. D. S. Naidu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 

and 


Dr. D. B. Price 
Spacecraft Controls Branch 
NASA Langley Research Center 
Hampton, Virginia, 23665 


Abstract: A common way of representing atmospheric density is by an 
exponential form using scale height, which is assumed to be constant over the 
whole interval of atmospheric altitude. In this simulation, the scale height 
has been readjusted depending upon the altitude Interval, and simulations are 
carried out for an aeroasslsted, coplanar orbital transfer vehicle. 


* See item 4 under List of Publications 
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TIME-SCALE SYNTHESIS OF A CLOSED- LOOP DISCRETE OPTIMAL CONTROL SYSTEM 


Dr. D. S. Naldu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA 23529 

and 


Dr. D. B. Price 
Spacecraft Controls Branch 
NASA Langley Research Center 
Hampton, VA 23665-5225 


Abstract: A two-time-scale discrete control system is considered. I he 

closed-loop optimal linear quadratic regulator for the system requires the 
solution of a full-order algebraic matrix Riccati equation. Alternatively, the 
original system is decomposed into reduced-order slow and fast subsystems. The 
closed-loop optimal control of the subsystems requires the solution of two 
algebraic Riccati equations of an order lower than that required for the 
full-order system. A composite, closed-loop suboptimal control Is created from 
the sum of the slow and fast feedback optimal controls. Numerical results 
obtained for an aircraft model show a very close agreement between the 
exact ( opt imal ) solutions and computationally simpler compos ite( subopt i mal ) 
solutions. The main advantage of the method is the considerable reduction in 
the overall computational requirements for the closed-loop optimal control of 
digital flight systems. 


* See item 6 under List of Publications 
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THREE-DIMENSIONAL ATMOSPHERIC ENTRY PROBLEM 

* 

USING METHOD OF MATCHED ASYMPTOTIC EXPANSIONS 


Dr. D. S. Naidu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 


Abstract: The analysis of a three-dimensional atmospheric entry problem using 
the method of matched asymptotic expansions is considered. A composite 
solution is formed In terms of an outer solution, an inner solution and a 
common solution. The outer solution is obtained from gravitationally dominant 
region, whereas the aerodynamical ly dominant region contributes to the inner 
solution. The common solution accounts for the overlap between the outer and 
inner regions. In comparison to the previous works, the present simplified 
methodology yields explicit analytical expressions for various components of 
the composite solution without resorting to any type of transcendental 
equations to be solved only by numerical methods. The method is applicable for 
obtaining autonomous guidance and control strategies for a variety of 
aerospace vehicles. 


* See items 9, 11 and 22 under List of Publications 
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FUEL-OPTIMAL TRAJECTORIES FOR AEROASSISTED 

* 

COPLANAR ORBITAL TRANSFER PROBLEM 


Dr. D. S. Naidu, Dr. J. L. Hlbey, and C. Charalambous 
Dept, of Electrical and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23529 


Abstract: We first describe briefly the various types of coplanar transfers. 
Then we address the fuel-optimal control problem arising in coplanar orbital 
transfer employing aeroassist technology. The maneuver involves a transfer 
from high Earth orbit to low Earth orbit and at the same time minimization of 
the fuel consumption for achieving the desired orbit transfer. It is known 
that a change in velocity, also called the characteristic velocity, is a 
convenient parameter to measure the fuel consumption. A suitable performance 
index is the total characteristic velocity which is the sum of the 
characteristic velocities for deorbit and for reorbit (or circularization). 

Use of Pontryagin minimum principle leads to a nonlinear, two-point boundary 
value problem in state and costate variables. The solution of the TPBVP is the 
stumbling block in obtaining fuel-optimal solution. This problem is solved by 
using a more efficient multiple shooting method which is a simultaneous 
application of a single shooting algorithm to equally divided points of the 
total interval of the solution. 


* See items 12, 17 and 25 under List of Publications 
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FUEL-OPTIMAL TRAJECTORIES OF AERO ASSISTED 

* 

ORBITAL TRANSFER WITH PLANE CHANGE 


Dr. D. S. Naidu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23529 


Abstract: The fuel-optimal control problem arising in noncoplanar orbital 
transfer employing aeroassist technology is addressed. The mission involves 
the transfer from high Earth orbit to low Earth orbit with plane change. The 
complete maneuver consists of a deorbit impulse to inject a vehicle from a 
circular orbit to elliptic orbit for the atmospheric entry, a boost impulse at 
the exit from the atmosphere for the vehicle to attain a desired orbital 
altitude and finally a reorbit Impulse to circularize the path of the vehicle. 
In order to minimize the total fuel consumption, a performance index is chosen 
as the sum of the deorbit, boost, and reorbit impulses. Application of 
optimization principles leads us to a nonlinear, two-point, boundary value 
problem, which is solved by using a multiple shooting method. 


* See items 13, 19, 20 and 29 under List of Publications 
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FUEL-OPTIMAL TRAJECTORIES FOR NONCOPLANAR 
ORBITAL TRANSFER VEHICLES 


Dr. D. S. Naidu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA 23529 


ABSTRACT: The fuel-optimal problem in noncoplanar orbital transfer employing 
aeroassist technology is addressed. The mission Involves the transfer from 
high Earth orbit to low Earth orbit with plane change. The complete maneuver 
consists of a deorbit impulse to inject a vehicle from a circular orbit to 
elliptic orbit to enter the atmosphere, a boost Impulse at the exit from the 
atmosphere for the vehicle to attain a desired orbital altitude and finally a 
reorbit impulse to circularize the path of the vehicle. In order to minimize 
the total fuel consumption, a performance index is chosen as the sum of* the 
deorbit, boost, and reorbl t impulses. For* a typical aeroassisted orbital 
transfer vehicle with high lift-to-drag ratio, the simulations are carried out 
using industry standard program to optimize simulated trajectories (POST). 


* See item 15 under List of Publications 
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SINGULAR PERTURBATIONS AND TIME SCALES 
IN THE DESIGN OF DIGITAL FLIGHT CONTROL SYSTEMS 


Dr. D. S. Naidu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA 23529 

and 

Dr. Douglas B. Price 
Spacecraft Controls Branch 
NASA Langley Research Center 
Hampton, VA 23665-5225 

Abstract: This paper Investigates the application of methodology of singular 
perturbations and time scales (SPaTS) to the control of digital flight 
systems. A block diagonal lzation method is developed to decouple a full-order, 
two-time (slow and fast) scale, discrete control system into reduced-order 
slow and fast subsystems. Basic properties and numerical aspects of the method 
are explored. This reveals an interesting fact that singularly perturbed 
discrete systems can be viewed as two-time scale systems. Next, the 
closed-loop optimal control of the two-time scale full-order system Involves 
the solution of a full order algebraic matrix Rlccati equation. Alternatively, 
using the block diagonal lzat ion method, the full-order system is decomposed 
into reduced-order slow and fast subsystems. The closed-loop optimal control 
of the subsystems requires the solution of only reduced-order algebraic matrix 
Rlccati equations. A composite closed-loop suboptimal control is constructed 
as the sum of the slow and fast optimal feedback controls. Numerical 
experimentation with an aircraft model shows close agreement between the exact 
solutions and the decoupled (or composite) solutions. The main advantage of 
the method is the considerable reduction in the overall onboard computational 
requirements for the evaluation of optimal guidance and control laws. It is 
believed that this paper also serves as a source of brief survey of digital 
flight systems. 


* See item 16 under List of Publications 
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NEIGHBORING OPTIMAL GUIDANCE FOR AN AEROASSISTED ORBITAL 
TRANSFER VEHICLE IN THE PRESENCE OF MODELLING UNCERTAINTIES 


Dr. D. S. Naidu, Dr. J. L. Hibey, and C. Charalambous 
Dept, of Electrical and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23529 


ABSTRACT: We intend to devise a neighboring optimal guidance scheme iur u. 
nonlinear dynamic system with stochastic Inputs and perfect measurements as 
applicable to fuel optimal control of an aeroasslsted orbital transfer 
vehicle. For the deterministic nonlinear dynamic system describing the 
atmospheric maneuver, a nominal trajectory is determined. Then, a neighboring, 
optimal guidance scheme is obtained, taking modelling uncertainties into 
account, a linear, stochastic, neighboring optimal guidance scheme is devised. 
Finally, the optimal trajectory is approximated as the sum of the 
deterministic nominal trajectory and the stochastic neighboring optimal 
solution. Numerical results are presented for a typical vehicle. 


* See items 23 and 27 under List of Publications 
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NEIGHBORING OPTIMAL GUIDANCE 

* 

FOR AEROASSISTED NONCOPLAN AR ORBITAL TRANSFER 


D. S. Naidu 

Dept, of Elect, and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 


Abstract: The fuel-optimal control problem in aeroassisted noncoplanar orbital 
transfer is addressed. The equations of motion for the atmospheric maneuver 
are nonlinear and the optimal (nominal) trajectory and control are obtained. 

In order to follow the nominal trajectory under actual conditions, a 
neighboring optimum guidance scheme is designed using linear quadratic 
regulator (LQR) theory for onboard real-time Implementation. One of the state 
variables is used as the independent variable in preference to the time. The 
weighting matrices in the performance index are chosen by a combination of a 
heuristic method and an optimal modal approach. The necessary feedback control 
law is obtained in order to minimize the deviations from the nominal 
conditions. 


* See item 24 under List of Publications 
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* 

ORBITAL PLANE CHANGE MANEUVER WITH AEROCRUISE 


Dr. D. S. Naldu 

Dept, of Electrical and Computer Engineering 
Old Dominion University 
Norfolk, VA 


Abstract: The synergistic plane change problem connected with orbital transfer 

employing aeroasslst technology, Is addressed. The mission involves transfer 
from high Earth orbit to low Earth orbit with plane change being performed 
within the atmosphere. The complete mission consists of a deorbit phase, 
atmospheric phase, and finally reorbit phase. The atmospheric maneuver is 
composed of an entry mode, a cruise mode, and finally an exit mode. During the 
cruise mode, constant altitude and velocity are maintained by means of bank 
angle control with constant thrust or thrust control with constant bank angle. 
Comparisons between these two control strategies bring out some interesting 
features. 


* See items 26 and 28 under List of Publications 
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PART II 


This final report under NASA Langley grant NAG1 736 consists of two 
parts. Part I consists of the summary of the earlier work whose reports have 
been submitted at various times during the period (1/1/87 to 12/31/1989). Part 
II describes in detail the research work done during the period 1/1/1990 to 


7/31/1990. 
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Abstract: The synergistic plane change problem connected with orbital transfer 
employing aeroassist technology, is addressed. The mission Involves transfer 
from high Earth orbit to low Earth orbit with plane change being performed 
within the atmosphere. The complete mission consists of a deorbit phase, 
atmospheric phase, and finally reorbit phase. The atmospheric maneuver is 
composed of an entry mode, a cruise mode, and finally an exit mode. During the 
cruise mode, constant altitude and velocity are maintained by means of bank 
angle control with constant thrust or thrust control with constant bank angle. 
Comparisons between these two control strategies bring out some interesting 
features. 


Nomenclature 


C : drag coefficient 

D 

C : drag coefficient at zero lift 

DO 

C : lift coefficient 

L 

C * lift coefficient for maximum lift-to-drag ratio 
LR 

D : drag force 

£ : maximum value of lift-to-drag ratio 

£ : aeropropulsi ve efficiency 

p 

g : gravitational acceleration 

H : altitude 

K : Hamiltonian 

j : specific fuel consumption 

sp 

i : inclination 

J : performance index 

K : induced drag factor 

L : lift force 

m : vehicle mass 

Q : Heating rate 

R : distance from Earth center to vehicle center of gravity 

R : radius of Earth 

E 

S : aerodynamic reference area 

T : thrust 

t : time 
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/ : velocity 

a : angle of attack 

/3 : inverse atmospheric scale height 

y : flight path angle 

5 : normalized density 

£ = normalized lift coefficient 

T) : thrust angle 

0 : down range angle 

X : costate (Langrange) variable 

\i : gravitational constant of Earth 

p : density 

<r : bank angle 

0 : cross range angle 

0 : heading angle 

AK : characteristic velocity 

Subscripts 

a : atmospheric boundary 

c : circularization at LEO 

d : deorbit at HEO 

e : entry to atmosphere 

f : exit from atmosphere 

j : beginning of aerocruise 

n : end of aerocruise 

s : surface level 
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I. INTRODUCTION 


The main function of space transportation system is to deliver payloads 
from Earth to various locations in space. Until now, this function has been 
performed by various rockets, the space shuttle, and expendable upper stages 
using solid or liquid propellants. In particular, considering the economic 
benefits and reusability, an orbital transfer vehicle (OTV) is proposed for 
transporting payloads between low Earth orbit (LEO) and high Earth orbit 
(HEO ). 1 The two basic operating modes contemplated for OTV are a ground-based 
OTV which returns to Earth after each mission and a space-based OTV which 
operates out of an orbiting hanger located at the proposed Space Station 
Freedom. 

In a typical mission, a space-based OTV, which is initially at the space 
station orbit (SSO), is required to transfer a payload to geosynchronous Earth 
orbit (GEO), pick up another payload, say a faulty satellite, and return to 
rendezvous with the orbiting hanger at SSO for refurbishment and redeployment 
of payload. The OTV on its return journey from GEO to SSO needs to dissipate 
some of its orbital energy. This can be accomplished by using an entirely 
propulsive (Hohmann) transfer in space only or a combination of propulsive 
transfer in space and aeroassisted maneuver in the atmosphere. It has been 
established that a significant fuel savings and hence increased payload 
capabilities can be achieved with propulsive and aeroassisted maneuvers 
instead of all-propulsive maneuvers . 1 This leads to an aeroassisted orbital 
transfer vehicle (AOTV), which on its return leg of the mission, dips into the 
Earth’s atmosphere, utilizes atmospheric drag to reduce orbital velocity and 
to achieve a desired orbital inclination. Basically, the AOTV performs a 
synergistic maneuver, employing a hybrid combination of propulsive maneuver in 
space and aerodynamic maneuver in the atmosphere. 

The plane change capability is required to (i) orbit a vehicle in a plane 
which does not pass through a launch site, (ii) shorten the time needed to 
reach multiple reconnaissance targets on a single orbital mission, (iii) 
reduce the time needed to return to base from orbit, (lv) perform effective 
rendezvous with satellites in different orbital planes, (v) avoid flights over 
hostile territory, and finally (vi) facilitate arrival and departure flights 
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from Space Station Freedom, in fulfilling specified mission objectives. 2 It 
should be noted that an orbital plane is usually defined in terms of 
inclination and longitude of the ascending node. For our present purpose only 
an inclination change is controlled. 

There are basically three methods of plane change, (1) impulsive method, 
(ii) aeroglide method, and (iii) aerocruise method. In impulsive method, the 
plane change is achieved entirely outside the atmosphere, and fuel consumption 
is prohibitively large for sizable changes of orbital plane. In both aeroglide 
and aerocruise methods, rockets are used to deflect the vehicle into the 
atmosphere, and the plane change is accomplished by heading change of the 
vehicle. With aeroglide there is no thrusting during the atmosphere, and with 
aerocruise, atmospheric drag is balanced by a continuous thrust to keep the 
spacecraft at a constant altitude and velocity. Propellent expenditure 
comparisons among the three methods of plane change show that the aerocruise 
method is superior to other competing methods for plane changes greater than 
about 20 degrees, and with heating restraints. The basic effect of propulsion 
during aerocruise is to (i) balance drag in order to maintain constant 
velocity, (ii) augment lift with a component of thrust, thus increasing 
cruising altitude over what it would be during aeroglide turn, and finally 
(iii) provide a lateral component of thrust giving the required turn necessary 
for plane change. The aeroglide and aerocruise methods utilizing atmospheric 

maneuver in conjunction with propulsion augmentation are also termed the 

3 4 

"synergistic” or "aeropropulsive" methods. 

, 5-7 

The following are some of the features of atmospheric plane change. 

(i) For plane changes of less than 15 degrees, an al 1 -propulsive maneuver is 
generally more efficient, (ii) An L/D of at least 2 Is required to offer a 
significant advantage over the all-propulsive plane change, and it is 
desirable to maximize the L/D of a vehicle, (iii) A plane change made at a 
node produces all inclination change whereas a turn at an orbit apex (90 
degrees from node) provides no inclination change, only a shift in the node. 
Hence, for maximum inclination change and minimum node shift, the turn should 
be centered over the node in the shortest possible duration. Thus, plane 
changes performed at maximum C l (l.e., high angle of attack) which are quicker 
are more fuel efficient than the slower maximum L/D turns, (iv) The total heat 
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load can be reduced substantially by carrying out a quicker high 
angle-of-attack turns rather than the slower maximum L/D turns, (v) An 
aerocruise (thrusting) turn offers significant advantages over an aeroglide 
(non-thrusting) turn, when the desired plane change is more than 10 degrees, 
(vi) During aerocruise, the high angle-of-attack and bank attitude of the 
vehicle produce a lateral component of thrust, which is responsible for a 
significant amount of plane change. 

Compared to other works, 5 ' 8 the highlights of the present work are (i) 

The analysis covers the complete mission from HEO to LEO. (11) The descent 
mode and ascent mode of atmospheric phase are analyzed with flight path angle 
as an independent variable, (iii) The time has been retained as an independent 
variable during the cruise mode, (iv) During the cruise mode, both bank angle 
control and thrust control are analyzed. 

In this report, we address the synergistic plane change problem arising 
in noncop lanar orbital transfer employing aeroassist technology. The mission 
involves the transfer from HEO to LEO with plane change being performed within 
the atmosphere. The complete mission consists of a deorbit phase, an 
atmospheric phase, and finally a reorbit phase. The atmospheric maneuver phase 
is composed of descent (entry) mode, cruise mode, and ascent (exit) mode. 
During the aerocruise mode, constant altitude and velocity are maintained 
either by ( i ) varying bank angle with constant thrust, or by (ii) varying 
thrust with constant bank angle. The comparison of these two control schemes 
bring out some interesting features. Numerical results are given for typical 

data. 
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II. MISSION DESCRIPTION 


For an orbital transfer problem, the following assumptions are made. (1) 
The initial HEO and final LEO orbits are circular, (ii) The mission is 
comprised of three impulses, (iii) The vehicle is represented as a constant 
point mass during atmospheric pass. Civ) A Newtonian inverse square 
gravitational field is used, (v) Earth's rotation is neglected, (vi) The 
atmosphere is exponential. 

The complete mission from HEO to LEO with atmospheric pass is depicted in 
Fig. 1. It consists of a deorbit phase, an atmospheric phase, and a reorbit 
phase. There are three impulses: first, a deorbit impulse at HEO to inject 

a vehicle into a HEO-entry elliptic orbit, second, a boost impulse AV b at the 
exit from the atmosphere for the vehicle to attain sufficient velocity to 
travel along an exit-LEO elliptic orbit, and finally, a circularizing Impulse 
AV to circularize the path of the vehicle. The atmospheric phase itself is 
composed of descent (entry) mode, cruise mode, and ascent (exit) mode. 

Consider the basic equations of motion for different phases of deorblt, 
aeroassist (or atmospheric flight), boost and reorbit (or circularization). 

Deorbit Phase 

Initially, we assume that a spacecraft is in a circular orbit of radius 
R , well outside the Earth’s atmosphere, moving with a circular velocity F^ 
Vp/R . Deorbit is performed by means of an impulse AF^, to transfer the 
vehicle from the circular orbit to elliptic orbit with perigee low enough to 
intersect the dense part of the atmosphere (Fig. ll. At D, since the elliptic 
velocity is less than the circular velocity, the impulse AF^ is executed so as 
to oppose the circular velocity F^. The deorbit Impulse AF^ causes the vehicle 
to enter the atmosphere of radius R ^ with a velocity F^ and flight path angle 
y It is known that the optimal-energy loss maneuver from the circular orbit 
is simply the Hohmann transfer and the impulse is parallel and opposite to the 
instantaneous velocity vector. 

Using the principle of conservation of energy and angular momentum at the 
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deorbit point D, and the atmospheric entry point E, we get, 
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V Z /2 - H/R = (V - A7 ) fz - f l/R 

c 1 a d d d 


(l) 


R V cos(-y ) = R (7 - A7 ) 

a e e da a 

from which solving for A7 , we get 

d 


A7 = -Jp/R 

d d 


2p( l/R - l/R )/[(R /R ) 2 /cos 2 y - l] 

-J a d ' L d a ' e J 


(3) 


It is easily seen that the minimum value of the deorbit impulse A7 d ^ 

obtained at y =0, corresponds to an ideal transfer wherein the space vehicle 
e 

grazes along the atmospheric boundary. To ensure proper atmospheric entry, the 
deorbit impulse A7^ must be higher than the minimum deorbit Impulse AT^ which 
is given by 


AF = A\i/R 

dm d 


2p( l/R - l/R ) /[( R /R ) 2 

^ a d /L da 


i] 


(4) 


Aeroassist (Atmospheric) Phase 

The atmospheric phase of the mission is composed of (i) descent mode, 
(ii) cruise mode, and (iii) ascent mode [Fig. 2]. 
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Ill DESCENT MODE 


During the descent mode, the equations of motion for the vehicle (without 


any thrusting) are given below [Fig. 3], The kinematic equations are, 
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3 ? = VBln * 


(5a) 


de 

dt 


= Vcos^cos\p/Rcos<p 


(5b) 


^ = Fcosrsini p/R 
at 


(5c) 


The force equations are 


dV 

ni-rr = - D - mgs l nr 
at 


(5d) 


mV^r = Lcosc + miV^R - g)cosy 
at 


(5e) 


mV ^ = i.sin<r/cosy - (m^^Ricosjrcosi/itan^ 

at 


(5f ] 


where, 


L = C (<x)pSV Z /2; D = C ( oc)pSV 2 /2 ; C = C + K(f 

L DU UU La 

g = I U/R 2 ; R s H + R ; p = p^expt-p//) 

Neglecting mass terms in comparison to aerodynamic terms in (5), we get 


dH 

dt 


= Fsiny 


(6a) 


dF 

dt 


^ = - D/m 


(6b) 


dy 

dt 


= Lcoscr/mF 


(6c) 
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^ = Lsln<r/mVcosy 
at 


(6d) 


Using flight path angle as the independent variable, we get 
dU 

?j- = 2mexp(/3H)siny/p SC coscr 
dr o l 


- C V/C COST 
dr d l 


(7b) 


-C- = tancr/cosr 

dr 


(7c) 


During the descent mode, let us assume that there is no banking of the 
vehicle, and hence no heading is achieved. Then (7) becomes 


— = 2mexp(/3H) si nr /p SC 
dr o l 


(8a) 


— = - C V/C 
dr d l 


Optimal Control Problem 


The optimal control problem is posed as follows. Given entry conditions, 
and the conditions at the end of the descent mode (or the Initiation of 
aerocruise mode), find the optimal control law which maximizes the final 
velocity, subject to altitude constraint H £ W fe . 11 This altitude constraint 
implies in a way heat-rate constraint. The performance index is given by 


The Hamiltonian for (8) and (9) is 


K = A 2mexp ( /3H) si nr /p SC 

H 0 L 


+ X 


- c v/c 

D 


L 


( 10 ) 


The adjoint equations are 
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(11a) 


dA /dy = -A 2ml3exp(l3H)siny/p SC 

H H 0 L 

dA /dr = A C /C 

v v D L 

Solving the state (8) and costate (11) equations, we get 
FA = V A = constant 

v e ve 


A exp(pH) = A exp(/3tf ) 

H ^ He e 


The boundary conditions for the adjoint variables are 


V* dV 


r = y, 


V r V a// 


r = y, 


= o 


From (12) and (13), 


A = 0 

H 


With (14), the Hamiltonian (10) reduces to 


K = - A V[C /C + K C ] 

v DO L L 


The optimal control is then given by 


o 

ac 

L 


leading to 


(lib) 


(12a) 

(12b) 


(13a) 


(13b) 


(14) 


(15) 


(16) 


C = <1 C ?K = C 

LO DO LE 


(17) 


where, C is the lift coefficient for maximum lift-to-drag ratio ((' /D ^ max = E 
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= 1/2-TkC . With optimal control (17) in (8), we solve for the 

DO 

altitude as 


F(y) = V exp[-(y-y )/£] 


H = In j^< 


exp (-0H ) + — 5 ^— (cosy - cosy ) 
r e p SC e 

0 LO 


1 - 1/0 


For small y, (19) reduces to 


H = ln|exp(-0/M - . ^ -(y 2 - y 2 ) 


e p SC 

K 0 LO 


-1/0 


At the start of the descent mode, y = y , and at the end of the 
(or the beginning of the aerocruise mode), y = = 0. Then the 


relations become 


"i ■ >"[■ 


ex p ( - pw e } + r§r 

0 LO 


-( 1 - cosy 


.»] 


- 1/0 


and with the approximation, 


H ^ = ln^exp(-0// e ) + - — 7 2 


p SC 

0 LO 


1/0 


Then, the inequality constraint on altitude H £ H , with (19a) 


y ^ cos 
e 


-1 


p SC 
* 1^1 1 


1 - - -P'-' 3 ".'} 


and with approximate solution (19b), 


y 2 5 

e 


p SC 

■ n 1 


r 1 

° - L — exp(-0tf ) - exp(-0// ) 

m(i [_ J « J 


velocity and 

(18) 

(19a) 


(19b) 

descent mode 
above 


(19c) 


( 19d) 

transforms to 
(20a) 


(20b) 
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The velocity at the end of the descent mode is obtained from (18) with T - 7 ^ 
= 0 as 


v } = ^ ex pC 2 v r V) 


( 21 ) 


13 


IV. AEROCRUI SE MODE: BANK ANGLE CONTROL 


We first write down the general equations of motion, inject the 
conditions for cruise flight, use the assumptions of small latitude, and 
finally optimize the heading change. During aerocruise mode, there is 
continuous thrusting. Thus the kinematic equations are [Fig. 3] 


dH 

dt 


Fsiny 


(22a) 


d0 

dt 


Vcosycos\p/Rcos<p 


(22b) 


^ = Vcosrsin \p/R 
dt 

The force equations are 


dV 

m-rr = T cost) - D - mgsinj 
at 


(22c) 


(22d) 


mV^r = (TsinT) + Ucoscr + m{V^R - g)cosy 
dt 


( 22e ) 


mV ^ = (Tsinr) + Usinir/cos? - ( mV?R) cos? cos\pta.n<p 
dt 


(22f ) 


The propulsion (thrusting) equation is 


% - - r/g/ 

at sp 


(22g) 


From the above equations of motion, we see clearly that during the 

atmospheric maneuver, if the lift vector L is rotated about the velocity 

vector V through the bank angle <r, it creates a lateral force component {Tsinri 

+ L)sincr orthogonal to the vertical plane that has the effect of changing the 

heading angle i//. At the end of the atmospheric phase, the equations (22c) and 

10 

(22f) for the cross range angle </>, and the heading angle tfr, become, 
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d<t>/dt _ _ tani p 
d\p/dt “ tan</> 


(23a) 


integration of which yields, 

cos<£cos0 = cosi (23b) 

where, i is the orbital inclination. For small values of cross range angle 0, 
the orbital inclination i is given by the heading angle 0 itself. Thus, the 
total change in the heading corresponds to the change in orbital inclination 
(plane change). 


Now let us insert the cruise conditions of constant altitude and 
velocity. 8 The constant altitude condition on (22a) gives zero flight path 
angle throughout. The constant velocity condition on (22d) boils down to 

Tcostj = D = pSV 2 C^{ a)/2 (24a) 

We note that the conditions at the beginning of the aerocruise mode are 
denoted by the subscript j. However, for simplicity in notation we shall 
continue to use the variables without any subscript to denote the cruise 
conditions. If the angle of attack a is held constant, then the drag force D 
is constant at a constant cruising altitude. Also, since the flight path angle 
is zero throughout, (22e) reduces to 

2 

(TsiriT) + £)coscr = m(g - V /ft) 

Combining (24a) and (24b), we get 

m(g - l^ 2 / ft) cost) 
cos<r - D[sinr) + (L/D)costj] 

where, 

slnr) + (L/D) cost) = (L + Tsim))/T = £ 


(24b) 


(25) 


(26) 
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Is called aeropropulsive efficiency. Alternatively, 


(Dtanr) + Dcoscr = m(g - V 2 /R) 


(27a) 


The above equation Is also rewritten as 


tamj = mK ^/ coscr - L/D 
where, = (g - V 2 /R)/D = ( gR - V )/RD 


(27b) 


(27c) 


From (24a), we see that for a given angle of attack a. If altitude H, and 
velocity V are kept constant, then the drag D and lift L forces are constant. 
The mass m always changes due to thrusting. Then the above relation (27b) for 
variable mass can be satisfied in any one of the following three ways. 


(i) Variable bank angle with constant angle of attack and thrust angle: With 
bank control, (22g) and (27a) mean that the thrust T Is constant leading to a 
constant mass flow rate. 

(ii) Variable thrust with constant bank angle: On the other hand, with thrust 
control, (24) implies that we need to change both magnitude and angle of the 
thrust, in order to keep a constant drag force. Thus for cruise condition, 
both thrust magnitude and angle need to be controlled such that Tcost) is 
constant, but Tsim) changes according to (24b) [see Fig. 3]. This leads to 
variable mass flow rate. 

(iii) Variable bank angle and variable thrust: Here, we change both bank angle 
and thrust magnitude and angle, in order to satisfy the cruise conditions 
(24a) and (24b). This also leads to variable mass flow rate. 


Obviously, the bank angle control leading to constant thrust (and hence 
constant mass flow rate) seems to be the simplest of all for implementation. 
However, it will be interesting to see which of the control schemes provides 
greatest amount of heading change and thereby inclination for the same amount 
of fuel expenditure. 
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The bank angle control with constant thrust magnitude and angle has been 
thoroughly discussed using arc length as independent variable. However, in 
our present work, we continue to use time as independent variable. Assuming 
the latitude to be small, the cruise motion is described by 


d0 

dt 


= Fcosi p/R 


< ^rr = ~ K /COST) 

dt 2 


= K tancr 
dt 3 


(28a) 


(28b) 


(28c) 


K = D/gl ; K = g/V - V/R < 28d) 

2 sp 3 

The bank angle control is given by 

coso* = mK /[tarn) + L/D] = mK (29a) 

1 4 

where, 

= (g - V 2 /R) _ 2 - 12 (29b) 

4 (L + Dtani)) pRSC o (a) [ tanr) + L/D] 

and k =VgR/V, the ratio of circular speed to cruise speed at R. From (29), we 
see that for a given angle of attack, and at constant altitude, speed, and 
thrust angle, the bank angle has to be varied as per the mass. That is, as 
mass 01 decreases along the flight, bank angle <r should be increased. Thus, in 
increasing the bank angle with the decrease of mass, we are trying to balance 
the decrease in the difference between the vehicle’s weight and centrifugal 
force with the sum of vertical components of lift and thrust. In this control 
scheme, both mass and bank angle change, whereas altitude, velocity, angle of 
attack, thrust, thrust angle, and mass flow rate are held constant. 


The cruise condition (29) reveals that 
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(a) With <r = 0°, there is no banking and the cruise conditions cam be 
maintained only by variable thrusting. 

(b) For 0 < cr < 90°, the cruising speed is less than circular speed. The lift 
is directed upward. The gravitational force is higher than the centrifugal 
force. 

(c) cr = 90° corresponds to cruise speed being equal to the circular speed, and 
all the aerodynamic force (TsinTj + L) is used for heading charge or turning. 
The gravitational force is equal to the centrifugal force. 

(d) For cr > 90°, the cruising speed is higher than the circular speed. The 
centrifugal force is higher than the gravitational force and hence the lift is 
directed downward in order to prevent the vehicle to escape from Earth. 

Given the initial values of mass and heading angle i{j , we find the 
initial bank angle <r from (29a). Also, we can solve (28b) directly as 

m(£) = - {K^/cosT))t + ( 30 ) 

Thus, cr( t ) , and \ pit) are solved until either of the desired final conditions 
m^, or i ft is realized. With a constant thrust angle 7), and given initial mass 
/n , and heading angle \p , and final mass m (or final heading angle \1) ) the 

J J n n 

sequence of solution of the aerocruise problem is to solve, first the mass 
equation (28b), second the bank angle equation (29a), and finally the heading 
angle equation (28c). 

In this formulation for aerocruise, we see that the heading angle changes 
with respect to bank angle as given by (28c), and bank angle in turn has to 
follow the mass as per (29a), and the mass varies independently according to 
(30). Hence, there is no optimization of heading angle w.r.t. bank angle 
control variable, for a given fuel consumption or of fuel consumption w.r.t. 
bank angle control variable, for a given heading change. Alternatively, from 
(29) and (30), we can solve for 0, which is now a function of thrust angle 77 , 
cruising altitude //, and cruising speed V . Then, we can find the optimal value 
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of 7 } , which should be maintained constant throughout the aerocruise to achieve 
maximum p. 

Optimization of Heading Change w.r.t. Thrust Angle 


The optimization problem here is to find an optimum thrust angle which is 
to kept constant throughout the cruise mode, in order to maximize the heading 
change. We can solve this problem in a variety of ways. Basically, the heading 
angle p can be solved from (28) to (30) in terms of time t, mass m, or bank 
angle <r. Thus, 


i/r = 0 + K InO /A ) + B 

n J 5|_ n J J 


B 

n 


where, 


J Sr 

K = 4Msi 

5 V L 


ini) + (L/D)cost) 


Ola) 


(31b) 


A = (1 + si no* )/coscr ; B = sincr ; 
l i i i l 


i = j, n 


(31c) 


Note that A , and B ^ can also be expressed in terms of mass m using (29), or 
in terms of time t with (30). For example, in terms of mass m t (31c) becomes 


[l + Jl - (m K ) 2 

L 1 4 J i 

A , = — J 7 ; B = h - [m K ) 2 ; i = J, n Old) 

i m a. l 14 

1 4 

In terms of bank angle cr, (31) is rewritten as, 11 


0 = 0, + 
n J 


I S 

sp 


sinr? + (L/D)costj 


] |H30s<r^(l + sincr K 

*[coscr (1 + sincr ) J 
^ n J J 


+ sincr - sincr 
J n 


(32) 


where, bank angle cr is related with thrust angle 7) as per (29). Now optimizing 

\p with respect to thrust angle tj, (i.e., making d\p /di) = 0) and assuming the 
n n 


initial 


value \1) = 0, 

J 


we get 
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cost} [ cost) + (L/0)sim)] 


(33) 


,-coscr (1 + sine ) 


ln< 


VI t SslIUT ;•» 

J "l 

cose (1 + sine ) f 

n j 


+ sine - sine 
J n 


= 0 


Using first order approximations in the change of the bank angle, 


cose = cos(e + Ae) = cose - (sine )Ae 
n J J J 


sine =sin(e + Ae) = sine + (cose )Ae 
n J J J 

and, linearizing the logarithm, the above transcendental equation 
becomes, 


(34a) 

(34b) 


(33) 


| cost) [cost) + (L/D)sinT)]sin 2 e^ + cos 2 e^ (1 + slne^)Ae = 0 


(35) 


Assuming that 1 + sine^ * 0, and Ae * 0, (35) becomes 


cost) [cost) + (L/D) sinT)]sin 2 e^ + cos 2 e^ = 0 


(36a) 


From (29a), 


cose = 


. - <. ; — 1—tk I A = m K 

j tarn) + L/D j l 


(36b) 


Using (36b) in (36a), 


[ tann + L/D] £( L/D) tan% + {(L/D) 2 - A 2 - l}tanT) - L/dJ 


(37) 


Again, assuming tarn) + L/D * 0, we finally get a simplified form as, 11 

(L/D) tan% + [(L/D) 2 - A 2 - l]tam) - L/D =0 (38) 

The implication of tarn) + L/D = 0 is that (a) 90° < T) < 180°, or (b) 270° < T) 
< 360°. The first condition implies that the thrust T will aid the drag force 
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instead of opposing it, and the second condition shows that the thrust T will 
oppose the lift instead of aiding it. 


Summarizing, for bank angle control we have 

^ = K tancr (39a) 

dt 3 


— « - K /cost) (39b) 

dt 2 

The cruise condition is given by 

coscr = mK^/[ tanTj + L/D] (39c) 

The sequence of solutions is 

(i) Using initial masses m ^ and L/D , solve (38) for optimal thrust angle tj. 

(ii) Using t?, fly and L/D, solve (39c) for ay 

(iii) Using tj, and m , solve for new m from (39b). 

(iv) Using the mass m , and tj, solve (39a) or (31) for new 0. 


(v) Go to (ii), and repeat the steps. 

(vi) Integration stops when m reaches the final m in (39b) and the 
corresponding maximum heading angle is obtained from (39a). 

n 

Although (39) can be solved for either set of given values of (a) m ^ and 
m or (b) </» and i b , the condition (38) requires m and m to determine the 

n J n J n 

optimal tj which is kept constant throughout the aerocruise mode. Hence, given 

the fuel consumption (m - m ) , we are trying to determine the maximum heading 

J n 

angle change. 
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V. AEROCRUISE MODE: THRUST CONTROL 


Here, we keep the bank angle constant throughout and change the thrust 

magnitude and angle in order to achieve the desired heading change and hence 
12 

the inclination. For the sake of simplicity, we repeat the equations at 
cruise, 


^ = Kcosi p/R 

(40a) 

dm „ , 

dt = ' K z /COSV 

(40b) 

= K tancr 
dt 3 

(40c) 

The cruise conditions are given by 


Tcost) = D 

(41a) 

(Tsinr? + Ocoscr = m{g - V 2 /R) 

(41b) 


Combining the two conditions, 

tanr? = mK^/cos<r - L/D (42) 

For a constant altitude //, speed F, and given angle of attack a and bank angle 
tr, as the mass m changes, the thrust angle r) follows (42). At the same time, 
the drag force D is to be kept constant as per (41a). Thus, in order to 
satisfy both the conditions (41), we need to adjust thrust magnitude T and 
angle rj in such a way that Tcost? is kept constant, and Tsinr) changes as per 
mass m. From (40), we see that the bank angle is kept constant throughout the 
cruise mode and hence the rate of change of heading angle is constant, whereas 
the mass flow rate is variable. This is in contrast to the bank angle control 
discussed in the last section. 
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Given a constant bank angle cr, the initial and final conditions , 

m (or 0 ), the sequence of solutions for the cruise flight with thrust 

n n 

control is first, solve the cruise condition equation (42) for t?, second solve 
mass rate equation (40b) for m , and finally solve the heading angle equation 
(40c) for 0. In the next section we try to find the optimal bank angle which 
should be kept constant throughout the cruise, to get maximum heading change. 


Optimization of Heading Angle w. r.t. Bank Angle 

Here, we are interested in finding the optimum bank angle so that the 
heading change is maximized. For this, we first solve (40)-(42) for the 
heading angle and then find the stationary value of ^ w.r.t. <r. Thus, 


0 

n 


(J g/F)sincrln 

sp 


secT) + 
J 

seer) + 

n 


tanrj 
_J 

tanT? 

n 


where, 


tanT? = m K /cos cr - L/D ; 
l i l 


2 

sec 7) 


= 1 + 


tan% 


i = J> 


n 


(43a) 


(43b) 


K = (7 g/V)sincr (43c) 

6 sp 

Alternatively, (43a) can be used to find the mass m for a given \f). Use of the 
stationary condition leads to 


cos crln 


’cost? ( 1 

+ sinT? )" 

J 

n 

COST) ( 1 

+ sinT) ) 

n 

J 


+ si 


2 (" 
no* 


(sin?) - si nr? ) + (L/D) (cost) - cost) 
n J n 


j ']- 0 


(44) 


Considering only the first order approximations in the thrust angle tj 

cost) = cos(t) + At)) = cost) - ( s i nT) )At) (45a) 

n J J J 

sinr) = sin(T) + At?) = sinT? + (cost? )At) (45b) 

n j J J 
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and linearizing the logarithm, the transcendental equation (44) becomes 


2 2 
cos cr + sin cr 


1 - ( L/D) tant? 


1 + tan T) 


\l 1 + sinr^ ] At? = 0 


Assuming 1 + sinr^ * 0; and At? * 0, (46) reduces to 


(46) 


2 2 
cos cr + sin cr 


1 - ( L/D) tanr? 


1 + tan 7) 


= 0 


(47) 


Using 


tarn? = j4/coscr - L/D 
J 


(48) 


in (47), we get a simplified quadratic equation in cos<r as 


( L/D) Acos^cr - [(L/D) 2 + 4 2 + 1 looser + {L/D) A = 0 


(49) 


It is interesting to note that the optimal conditions (38) and (49) for bank 
control and thrust control respectively, are interchangeable by the cruise 
condition (36b) or (48), 


Summarizing, for thrust control we have 
^ = K tan o- (50a) 

dt 3 

= ~ K /cost) (50b) 

dt 2 

The cruise conditions are given by 


7 COST) = D 

(50c) 

tan^ = mK^/coscr - L/D 

(50d) 
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The sequence of solutions is 


(i) Using initial mass m ^ and L/D, solve (49) for optimal bank angle < r. 

(ii) Using c r, and L/D, solve (50d) for tj^. 

(iii) Using T) ^ , find the thrust T ^ from (50c). 

(iv) Using i /r , and <r, solve (50) for new 0. 

(v) Using 17 ^ , and m j( solve for new m from (50b) or (43). 

(vi) Go to (ii), and repeat the steps. 


(vii) Integration stops when m reaches the final m in (50b) and the 
corresponding maximum heading angle 0^ is obtained from (50a). 

Although (50) can be solved for either set of given values of (a) m and 
m or (b) i/i and 0 , the condition (49) requires m and m to determine the 

n J n J n 

optimal o- which is kept constant throughout the aerocruise mode. Hence, given 

the fuel consumption (in — tn ) , we are trying to determine the maximum heading 

J n 

angle change. 
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VI. ASCENT MODE 


This is just a replica of the descent mode except for the change in the 
mass of the vehicle, and the boundary conditions. Thus, we go through the 
equations of motion, change the independent variable to flight path angle, and 
finally assume that there is no appreciable change in heading angle. Thus, 


^ = 2m exp (pH) si nr /p SC 
dy n 0 L 


(51a) 


= - C V/C (51b) 

dr d l 

The optimal control problem is posed as follows. Given initial conditions (or 
the conditions at the end of the cruise mode), and the conditions at the end 
of the ascent mode, find the optimal control law which maximizes the final 
velocity. 10 The performance index is given by 

J = - V (52) 

f 

The Hamiltonian for (51) and (52) is 


H = A 2m exp (0//) si nr /p SC 

H n 0 L 


+ A 


- c V/C 

D L_ 


The adjoint equations are 


dX /dr = -A 2m /3exp(/3//)sinr/p SC 

H H n 0 L 


(53) 


(54a) 


dA /dr = A C /C (54b) 

v y D L 

Solving the state (51) and costate (54) equations, we get 

l^A = V A = constant (55a) 

v n vn 


26 


A exp(/3h) = A exp(/3// ) 

H Hn n 


(55b) 


The boundary conditions for the adjoint variables are 


\ r , dJ 

A (y = y) = oTT 
v r oK 


y = r. 


= - 1 


(56a) 


\ f , dJ 

V r = = 3W 


r = y 


= o 


(56b) 


From (55) and (56), 


A = 0 

H 


(57) 


With (57), the Hamiltonian (53) reduces to 


K = - A V[C /C + K C ) 

v DO L L 


(58) 


The optimal control is then given by 


™ = 0 

dC 

L 


(59) 


leading to 


C = 4c /K = c 

LO DO LE 


(60) 


where, C is the lift coefficient for maximum lift-to-drag ratio [L/D ) = E 

LE max 

= 1/2 TKc . With optimal control (60) in (51), and noting that y = 0, we 

DO n 

solve for the velocity and altitude as 


Wy) = V exp (-y/£) 

n 


(61) 


H = In 


exp(-/3H ) 

n 


2/3/7? 

^sc :L<cos, ' 

f r r\ 



(62a) 
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For small y, (62a) reduces to 


H = In 


(3m 

exp(-(3H ) - c ~ n - y 2 

n P oO 

0 LO 


- 1/0 


(62b) 


At the end of the ascent mode, y = y and H 


H . Then (62) becomes 
r 


H = In 

f 


0m 

. \ n 2 

exp( —(3H ) - j~sc-7 f 

0 LO 


- 1/0 


(63) 


Using (19d) and (63) and noting that = H f and ffj = H n> we get 


y = - y m /m 

f eN J n 

The relation (64) shows that at the end of the atmospheric phase, the vehicle 

has to leave the atmosphere with a positive flight path angle higher in 

magnitude to that of the entry flight path angle. This is due to the fact that 

the mass m at the beginning of the ascent mode (or end of the cruise mode) is 
n 

less than the mass m^ at the end of the descent mode (or beginning of the 
cruise mode). 
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VII. BOOST AND REORBIT PHASE 


During the atmospheric flight, the vehicle performs the desired plane 

change and dissipates some energy due to atmospheric drag. Therefore, a second 

Impulse is required to boost the vehicle back to orbital altitude. The vehicle 

exits the atmosphere at point F, with a velocity and flight path angle 7 f • 

The additional impulse AF b , required at the exit point F for boosting into an 

elliptic orbit with apogee radius R_ and the reorbit Impulse required to 

insert the vehicle into a circular orbit at point C, are obtained by using the 

principle of conservation of energy and angular momentum at the exit point F, 

n 

and the circularization point C. Thus, we have, 

( V + AV ) Z /2 - n/R = [V - tV ) 2 /2 " M/R (65) 

f b ' a c c c 

{V + A V )R cos^ = R {V - Atf ) 

f b a f c c c 

Solving for AV and AV from the above equations (65) and (66), 
b c 


AV = 1 2p( 1/R - 1/R )/[l ~ (R /R ) 2 cos 2 y ] - V (67) 

b^ a c' 1 - ac t i 

AV = ^7r" - |2p(l/R - 1/R )/[(R /R ) 2 / cos \ ■ ^ (68) 

c c i a c c a 

Finally, the vehicle is in a circular orbit (of radius R c ) moving with 

the velocity V = V\i / R , 
c c 
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VIII. NUMERICAL DATA AND RESULTS 


The following set of data is used for a typical orbital maneuvering 
research vehicle. 

Orbital Data 

Altitude of HEO, H = 115,000 m 

d 

Altitude of LEO, H = 115,000 m 

c 

Altitude of atmospheric boundary, H* - 110,000 m 
Radius of Earth, R = 6,356,766 m 

E 

2 

Acceleration due to gravity at sea level, g Q = 9.80665 m/sec 

3 

Atmospheric density at sea level, p Q = 1.225 kg/m 
Gravitational constant of Earth, p = 3.986x10 m /sec 
Inverse atmospheric scale height, (3 = 1/7280 m 

Vehicle Data 

Initial mass, m = 4760 kg 

Propellent available for cruise, (m - m ) = 1810 kg 
K J n 

Final mass, m = 2950 kg 

n 

2 

Aerodynamic reference area, S = 11.613 m 
Specific fuel impulse, I ^ = 290 sec 

The aerodynamic characteristics are described in terms of the angle of attack 
a as 

C = -2.068686996a 3 + 2.943200144a 2 + 0.080347684a + 0.031320026 (69a) 

L 

C = 0.267339707a 3 + 1.814473159a 2 - 0.389985867a + 0.068372034 (69b) 

D 
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Fig. 4 shows the variations of the lift and drag coefficients C^, and C D and 
Fig. 5 shows the variations of the lift-to-drag ratio E, and the 
aeropropulsive efficiency E as a function of the angle of attack a. The 
maximum lift-to-drag ratio of 2.3149 occurs at the angle of attack of 13 
degrees. 

Deorbit Phase 

Initially, the vehicle is at a HEO altitude of 115 km orbiting with a 

circular velocity V of 7847.97 m/sec. A deorbit impulse LV ^ of 518.99 m/sec 

puts the vehicle in an elliptic orbit to intersect the atmospheric boundary at 

an altitude H of 110 km. At the atmospheric entry point, the velocity V is 
a 

7334.17 m/sec and the flight path angle r e Is -0.77 deg. 

Atmospheric Phase: Deorbit Mode 

During the descent mode of atmospheric phase, the vehicle descends from 
an altitude H of 110,000 m at a velocity V of 7334.17 m/sec to a cruise 

a * 

altitude H of 72,521 m and the cruise velocity V of 7291.17 m/sec according 
to relations (18) and (19). During this time, the vehicle is maintained at an 
angle of attack of 13 degrees corresponding to the lift given by (17) for 
maximum lift-to-drag ratio. The time solutions of altitude, velocity, and 
flight path angle for the descent mode are shown in Fig. 6. The time taken for 
descent mode is found to be 494 seconds. 

Cruise Mode: Bank Angle Control 


The cruise mode is analyzed using bank angle control or thrust control. 
With bank angle control, for the same fuel consumption and a given L/D (or 
angle of attack a), the optimum thrust angle as obtained from (38) is 
represented in Fig. 7. With this optimum thrust angle, the corresponding 
heading angle is obtained from (32) and is shown in Fig. 8. It is seen that at 
these constant cruising conditions (of altitude of 72521 m, velocity of 7291.7 
m/sec, and thrust of 6108.9 Nw) , the maximum heading and hence maximum 
inclination of 18.56 degrees is achieved with a thrust angle of 44.72 degrees 
and at a higher angle of attack of 24 degrees rather than at the angle of 
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5“8 

attack of 13 degrees corresponding to maximum L/D. Corresponding to the 
maximum inclination of 18.56 degrees, the time solutions for mass, bank angle, 
and heading angle are shown in Fig. 9, where the total time taken for the 
cruise mode is 830 seconds. 

Cruise Mode: Thrust Control 

With thrust control, for the same fuel consumption and a given L/D (or 
angle of attack a), the optimum bank angle as obtained from (49) is shown in 
Fig. 10. With this optimum bank angle, the corresponding heading angle is 
obtained from (43) and is shown in Fig. 11. It is seen that at these cruising 
conditions, the maximum heading angle of 17.7 degrees is achieved with a bank 
angle of 51.9 degrees and at a higher angle of attack of 20 degrees rather 
than at the angle of attack of 13 degrees corresponding to maximum L/D. 
Corresponding to the maximum heading angle of 17.7 degrees, the time solutions 
for mass, thrust, thrust angle, and heading angle are shown in Fig. 12, where 
the total time taken for cruise mode is 1345 seconds. 

The comparison of maximum heading angle as a function of angle of attack 
for both the control strategies shown in Fig. 13, indicates the superiority 
of bank control over thrust control. It is to be noted that the heading 
achieved depends on the type of control used, cruise conditions, and the angle 
of attack. 

Atmospheric Phase: Ascent Mode 

At the end of cruise mode, the vehicle ascends to the atmospheric 
boundary with a constant angle of attack of 13 degrees corresponding to 
maximum lift-to-drag ratio as given by (60). At the end of the ascent mode, 
the exit velocity is 7238.1 m/sec, the flight path angle as given by 
(64) is 0.9781 deg. The time solutions are shown in Fig. 14. 

For the atmospheric phase with bank angle control, the total solutions 
for altitude, velocity, flight path angle, heading angle, and heating rate are 
shown in Fig. 15. Similarly, total solutions for thrust control are shown In 

5 

Fig. 16. The heating rate is computed from 
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Q = 3.08xlO'V / V 3,08 Watts/cm 2 (70) 

* k k 

3 

where , p is expressed in kg/km , and V is expressed in km/sec. 

k k 

Reorbit Phase 

At the end of the atmospheric phase, a boost impulse LV of 380 m/sec is 

b 

executed to bring the vehicle to its original altitude H of 110 km. At this 
time, once again a circularizing impulse Atf of 247.97 m/sec is imparted to 

c 

finally put the vehicle in circular orbit. 
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IX. CONCLUDING REMARKS 


We have addressed the synergistic plane change problem in connection with 
orbital transfer employing aeroassist technology. The mission Involved 
transfer from high Earth orbit to low Earth orbit with plane change being 
performed within the atmosphere. The complete mission consisted of a deorbit 
phase, an atmospheric phase, and finally a reorbit phase. The atmospheric 
maneuver was composed of an entry mode, a cruise mode, and finally an exit 
mode. The descent and ascent modes have been analyzed using flight path angle 
as an independent variable for maximizing the cruise and exit velocities with 
a constraint on the minimum cruise altitude. During the cruise mode, constant 
altitude and velocity were maintained by means of bank angle control with 
constant thrust or thrust control with constant bank angle. Conditions have 
been obtained for maximizing the heading angle. Under given cruising 
conditions, the maximum heading angle has been achieved with an angle of 
attack higher than that corresponding to the maximum llft-to-drag ratio. 
Comparison between the two control strategies has shown the superiority of 
bank control over thrust control in terms of the maximum achievable heading 
angle. 
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Fig. 3. Coordinate System 
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Fig. 5. Lift-to-Drag Ratio and Aeropropulsi ve Efficiency 
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Fig. 6(a). Altitude as a Function of Time during Descent Mode 
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Tig, 6(b). Velocity as a Function of Time during Descent Mode 
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Fig. 6(c). Flight Path Angle as a Function of Time During Descent Mode 
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Fig. 8. Maximum He 
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Fig. 9(a). Mass as a Function of Time with Bank Angle Control 
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Fig. 9(b). Bank Angle as a Function of Time with Bank Angle Control 
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Fig. 12(a) 
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ig. 12(b). Thrust as a Function of Time with Thrust Contr 
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ATMOSPHERIC PHASE WITH BONK CONTROL 
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Fig. 15(b). Velocity as a Function of Tine with Bank Angle Control 
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Fig. 15(c). Flight Path Angle as a Function of Time with Bank Angle Control 


0000000000 



63 


Fig. 15 (d). Heading Angle as a Function of Time with Bank Angle Control 
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Fig. 15(e). Heat Rate as a Function of Time with Bank Angle Control 
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Function of Time with Thrust Control 
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Fig. 
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Fig. 16(d). Heading Angle as a Function of Time with Thrust Control 
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Fig. 16(e). Heat Rate as a Function of Time with Thrust Control 


